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ABSTRACT 

MHD in protostellar discs is modified by the Hall current when the ambipolar diffusion 
approximation breaks down. Here I examine the Balbus-Hawley (magnetorotational) 
instability of a weak, vertical magnetic field within a weakly-ionized disc. Vertical 
stratification is neglected, and a linear analysis is undertaken for the case that the 
wave vector of the perturbation is parallel to the magnetic field. 

The growth rate depends on whether the initial magnetic field is parallel or an- 
tiparallel to the angular momentum of the disc. The parallel case is less (more) unsta- 
ble than the antiparallel case if the Hall current is dominated by negative (positive) 
species. The less-unstable orientation is stable for \ ^ 0.5, where \ is the ratio of a 
generalised neutral-ion collision frequency to the Keplerian frequency. The other ori- 
entation has a formal growth rate of order the Keplerian angular frequency even in the 
limit x - * 0! In this limit the wavelength of the fastest growing mode tends to infinity, 
so the minimum level of ionization for instability is determined by the requirement 
that a wavelength fit within a disc scale height. In the ambipolar diffusion case, this 
requires x > v a/cs] in the Hall case this imposes a potentially much weaker limit, 
X > v 2 A /c 2 s . 
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1 INTRODUCTION 



The Balbus-Hawley (magnetorotational) instability (Ve- 
likhov 1959; Chandrasekhar 1961) of weak magnetic fields in 
accretion discs drives MHD turbulence that transports an- 
gular momentum radially outwards (Balbus & Hawley 1991, 
Hawley & Balbus 1991; Stone et al 1996), and is thought to 
play an important role in the evolution and dynamics of as- 
trophysical accretion discs. The instability has also been in- 
voked as a component of a disc dynamo model, in which the 
instability creates radial field from vertical field, the shear in 
the disc creates azimuthal field from the radial component, 
and the Parker instability creates vertical from azimuthal 
field and expels flux from the disc (Tout & Pringle 1992) . 

If the disc is sufficiently ionized that flux-freezing is a 
good approximation, the characteristic growth rate and ver- 
tical wave number of the instability are of order SI and SI/va 
respectively, where Q and va are the Keplerian angular fre- 
quency and the Alfven speed in the disc. The growth of the 
instability in a weakly-ionized, magnetized disc is of interest 
for theories of star formation and the subsequent evolution 
of protostellar discs. In this context, the degree of coupling 
between the field and neutral gas is determined by the trace 
charged species that are produced by cosmic-ray, radioac- 
tive or thermal ionization (Hayashi 1981) or the X-ray flux 
from the central, magnetically active star (Glassgold, Najita 
& Igea 1997). The level of ionization may only be sufficient 



to couple a magnetic field to the material in the surface re- 
gions, where incident X-rays and cosmic rays are absorbed, 
restricting magnetic activity to the surface layers over a large 
portion of the disc (Gammie 1996; Wardle 1997). 

Two models for field diffusion are generally adopted 
when considering the breakdown of flux freezing in weakly- 
ionized gas. At relatively low densities, the magnetic field 
can be regarded as frozen into the ionized component of the 
gas and drifts with it through the neutrals, a process re- 
ferred to as ambipolar diffusion (Spitzer 1978). The linear 
growth of the instability has been examined in the pres- 
ence of ambipolar diffusion (Blaes & Balbus 1994), and the 
nonlinear development has been investigated (MacLow et al 
1995; Hawley & Stone 1998). The instability grows only if 
the neutral- ion coupling time scale is less than about Q -1 . 
At high densities collisions with neutral particles stop the 
charged species drifting with the magnetic field, and the 
fluid can be treated as resistive - this is the limit of Ohmic 
diffusion. The behaviour of the instability in the presence 
of Ohmic diffusion has also been considered in the linear 
(Jin 1996) and nonlinear regimes (Sano, Inutsuka & Miyama 
1998). In this case it is the resistive diffusion time scale that 
has to be compared with fi . 

However, both ambipolar and resistive diffusion are gen- 
erally poor approximations in protostellar discs, where the 
gas density is sufficient for collisions with neutrals to par- 
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tially or completely decouple ions and grains from the mag- 
netic field, but electrons are either well-coupled or partially 
coupled to the field (Wardle & Konigl 1993). Instead, the 
gas is in the Hall regime, where variations in the decoupling 
amongst the charged species produces an overall handedness 
in the fluid with respect to the magnetic field (Wardle & Ng 
1999). That is, the fluid dynamics is no longer invariant un- 
der a global reversal of the magnetic field. For example, left- 
and right-circularly polarized Alfven waves propagate at dif- 
ferent speeds, and damp at different rates (Pilipp, Hartquist 
& Havnes 1987; Wardle & Ng 1999). This qualitative change 
in the wave modes supported by the fluid implies that the 
dynamical evolution of the magnetized gas will be dramati- 
cally modified. 

This paper examines the modifications to the Balbus- 
Hawley instability in weakly-ionized gas. In Section 2 I write 
down the fluid equations appropriate for a weakly-ionized, 
magnetized, near-Keplerian disc. The linearized equations 
and dispersion relation are obtained in Section 3 for pertu- 
bations with wave vector parallel to an initially vertical field, 
neglecting stratification of the disc. Two important dimen- 
sionlcss parameters emerge from this analysis: x> the ratio of 
the neutral coupling frequency to the Keplerian frequency, 
and the ratio of the Hall and Pedersen conductivities, ai and 
<72- The dependence of the growth rate on these parameters 
is presented in Section 4, where I show that the growth rate 
of the instability in the Hall regime is dependent on whether 
the initial vertical field is parallel or anti-parallel to the disc 
angular velocity vector. In particular, the surprising result 
emerges that the growth rate may be non-zero even in the 
limit x ~~ * 0. This result is interpreted physically using a 
two-fluid model of the Hall limit in Section 5. The impli- 
cations for the evolution and dynamics of protostellar discs 
are discussed in Section 6, and the paper is summarized in 
Section 7. 



2 FORMULATION 

I begin by writing down the equations describing the 
fluid dynamics of a weakly-ionized, magnetized, non-self- 
gravitating disc moving in the potential of a central point 
mass M. The definition of weakly-ionized in this context 
is that the inertia of the charged species can be neglected. 
This constrains the frequencies of interest to be below the 
collision frequency of each charged species with the neutrals. 
Under this assumption the charged species can be incorpo- 
rated into a conductivity tensor a , and seperate equations 
of motion are not needed. 

Following standard practice, I write the equations in 
a local Keplerian frame, so v is the departure from Kep- 
lerian motion - the fluid velocity in the laboratory frame 
is v + vk , where the Keplerian velocity vk ~ \J GM/r <j> 
in the canonical cylindrical coordinate system (r,cf>,z). The 
partial time derivatives in the Keplerian frame, written as 
d/dt in the equations to follow, correspond to + Q-§^ in 
the laboratory frame. The continuity equation is 



| + v.M = o. 



(i) 



is cancelled by the portion of the centripetal term associated 
with exact Keplerian motion. Tidal effects can be neglected 
and the momentum equation becomes 

^ + (v ■ V)« + -VP - 2Q.V4,t + ±nv r 4> = 12LJL ( 2 ) 
at p z cp 

where Q = u^/r is the Keplerian frequency. The pressure 
gradient term does not appear in the linearized equations, 
which are restricted to Alfven modes, thus we need not con- 
sider an equation of state. The current density is given by 
Ampere's law 

J=-^VxB. (3) 

47T 

The evolution of the magnetic field B is determined by the 
induction equation 



dB 

at 



= V x (v x B) - cV x E' - lQ.Br 



(4) 



where E' is the electric field in the frame comoving with the 
neutrals, and B must satisfy 



V- B = 0. 



(•») 



The last term in eq. (^) represents the generation of toroidal 
field from the radial component by Keplerian shear. 

In the weakly-ionized limit relevant to protostellar discs, 
J and E' are related by a conductivity tensor, a, that is de- 
termined by the abundances and drifts of the charged species 
through the neutral gas. I characterize each charged species 
j by particle mass rrij and charge Zje, number density rij, 
and drift velocity through the neutral gas Vj, and assume 
overall charge neutrality, V^.rijZj = 0. Assume that the 
fluid evolves on a time scale that is long compared to the 
collision time scale of any type of charged particle with the 
neutrals. Then each charged particle drifts through the neu- 
trals at a rate and direction determined by the instantaneous 
Lorentz force on the particle and the drag force contributed 
by neutral collisions: 



(e' 



x B 



- IjmjpVj = , 



(6) 



where —< av >j /(m + rrij), < av >j is the rate coeffi- 
cient for momentum transfer by collisions with the neutrals 
and m is the mean neutral particle mass. The relative impor- 
tance of the Lorentz and drag forces in balancing the electric 
force on species j is determined by the Hall parameter, 



Pj — 

rrijC yjp 



(7) 



Near the disc midplane, on scales small compared to the disc 
thickness, the radial gradient in the gravitational potential 



which is the product of the gyrofrequency and the time 
scale for momentum exchange by collisions with the neu- 
trals. This parameter controls the magnitude and direction 
of the charged particle drift and whether a particular species 
may be regarded as tied to the magnetic field by electromag- 
netic stresses (\/3j \ 3> 1 implies cE' w — Vj X B), to the neu- 
tral gas by collisions (\(3j\ <C 1 implies "fjinjpVj w ZjeE'), 
or to neither (\/3j\ ~ 1). 

The current density can be expressed in terms of E' and 
B by inverting eq ([]) for Vj, and forming J = rijZjVj. 
The resultant expression can be written as 

J — <T||23jj + (j\B x E'j_ + <72-Ej_ (8) 

where E[, and E'j_ are the decomposition of E' into vectors 
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parallel and perpendicular to B respectively and the com- 
ponents of <r are the conductivity parallel to the magnetic 
field, 



rijZ 2 e 2 ec 



mjjjp B 
i 3 

the Hall conductivity, 



ec x — *s 

= g } , njZjfy 



B 



72 7 



and the Pedersen conductivity 

1 ec \ ^ n,jZj(3j 

1 + 0? = 1 + 



(7_> 



En,i7 2 e 2 



(9) 



(10) 



(11) 



(Cowling 1957; Norman & Heyvaerts 1985; Nakano & Ume- 
bayashi 1986). Later I shall refer to the total conductivity 
perpendicular to the field, 



cr_L = \/ o\ + a 



(12) 



The second forms given for the parallel and Pedersen 
components in eqs (^) and ( pi] ) are useful for comparing 
their magnitudes with the Hall component, cr n and 02 are 
always positive, whereas o\ may take on either sign, de- 
pending how the magnitudes of the Hall parameters of dif- 
ferent charged species are distributed. The ambipolar dif- 
fusion limit is recovered when |/3j| 2> 1 for most species, 
when Cm 3> 02 2> |<xi|. The Ohmic limit is approached when 
\j3j\ -C 1 for most species, implying that 02 ~ <7|| 2> |ci|. 
Less well-known is the Hall limit, which can be approached 
when the magnitudes of the Hall parameters associated with 
particles of one sign of charge are large, and those of the 
other sign are small. Then |cti| is greater than both cry and 
(J2 (see Section]^). 

Numerical evaluation of the Hall parameters for ions 
and electrons (e.g. Wardle & Ng 1999) gives /3; ~ 
0.46B 

mG^iQ and (3 & ~ — 3507^q -B m G^\oo^ > where B m Q is 
the magnetic field strength in mG, nio and T100 are the den- 
sity of hydrogen nuclei and electron temperature normalised 
by 10 10 cm -3 and 100 K respectively. For all but the smallest 
grains, one finds \/3 g \ <C 1, on account of their large geomet- 
ric cross-section. The role played by grains is uncertain as 
they rapidly settle to the disc midplane. As the gas den- 
sity at a particular cylindrical radius varies by many orders 
of magnitude from the disc midplane to the surface, and 
the magnetic field strength may range from the milliGauss 
strengths found in the parent cloud to values a hundred 
times greater, one expects large regions of protostellar discs 
to be in a regime intermediate between the ambipolar diffu- 
sion and resistive limits, where o\ is non-negligible (Wardle 
& Ng 1999). 



3 LINEARIZATION 

The equations are linearized about an initial steady state in 
which v, VP, E' and J vanish, and B — Bz. This is ap- 
propriate for a box that is small in all dimensions compared 
to the scale height. Changes in the conductivity associated 
with the perturbations do not appear in the linearized equa- 
tions, as E' vanishes in the unperturbed state. Thus I need 



not specify how the charged particle abundances respond to 
the perturbations - all that is needed are the components of 
the conductivity tensor in the unperturbed state. 

Assuming solutions to the linearized equations of the 
form expi(ut — kz), the linear system splits into one sub- 
system corresponding to sound waves propagating along the 
magnetic field, and the interesting subsystem, for which the 
^-components of the perturbed velocity, magnetic field, cur- 
rent density and electric field all vanish. Denoting the result- 
ing two-dimensional vectors by the subscript _L, one finds the 
following expressions for the perturbations. 

Ampere's law allows SJ± to be expressed as 



c T ick 
SJ± = — - 

4n 



-1 




SB 



(13) 



where the matrix on the right hand side represents the cross 
product operator z X . Substituting this into the linearized 
momentum equation allows the perturbations in v and B 
to be related: 



2iQ 



Sv ± , SB 

= KVA 



1>A 



B 



(14) 



In the absence of rotation (Q — 0), this reduces to 
the standard relationship for Alfven waves, SB±/B — 
(to /kvA)Sv j_ / 'va- The ratio to/kvA is 1 in the ideal case, 
but is modified in the non-ideal case (Wardle & Ng 1999). 
The off-diagonal terms in the matrix on the left hand side 
are introduced by the Coriolis and centripetal terms in eq. 
(^), implying that the unstable mode is a modified Alfven 
mode. 

The linearized induction equation becomes, after sub- 
stituting for Sv± , 



uj 2 - k 2 v\ + 3fi 2 
— 2if2o; 



2iQ,u> 
j 2 - k 2 v\ 



SB, 



ike 



29. 



8E\ 



(15) 



The effects of finite conductivity are contained on the right 
hand side. In ideal MHD, SE' ± — 0, and the dispersion 
relation is obtained by setting the determinant of the matrix 
on the left to zero. In the absence of rotation, this leads to 
the dispersion relation for Alfven waves; with SI 7^ the 
dispersion relation for the magnetorotational instability is 
obtained (cf. Balbus & Hawley 1991). 

The conductivity tensor relates SE' ± to SJ± , and hence 
via ( |l^ ) to another relationship between SB± and SE' ± : 



SE\ =- 



where 



ike 
4na 2 



Sai —(72 
(72 S(7l 



SB 



s = sign(B z ). 



(16) 



(17) 



Ideal MHD is a good approximation for long wavelength 
perturbations, when ( |16[ ) shows that \SE\ \ is small. It 
breaks down when \SE' ± \ implied by ( |16[ ) becomes signif- 
icant in the linearized induction equation ([Hs]) , that is when 
k 2 c 2 /(Ana±) <; |w|. This corresponds to kvA ii i^c, where 
the critical frequency 



B 2 (7J 

pc 2 



(18) 
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is the generalisation of the ion-neutral collision frequency in 
the ambipolar diffusion case (Wardle & Ng 1999). 

The strength of the coupling between the magnetic field 
and the disc is determined by the parameter 

*=7P ^ 

which is much larger or less than unity if the field is well 
or poorly coupled to the disc respectively. \ reduces to the 
parameter 'ypi/il that appears in the ambipolar diffusion 
case (Wardle & Konigl 1993; Blaes & Balbus 1994). If the 
diffusion is Ohmic, then x = ^/fi c , where Q, c — k 2 c 2 / A-ko 2 
is the diffusion rate for waves with k = Q/va (cf. Jin 1996). 

Substituting for 8E' ± in ( |l5| ) yields a matrix whose de- 
terminant must vanish if nonzero solutions for the pertur- 
bations are to exist. The resulting dispersion relation is of 
the form 



a{kv A /Q) A + b(kv A /Q) 2 + c = 



(20) 



where 

a = X 2 + (lso 1 +2o 2 v)cTl 1 x + l + v 2 , (21) 

b = (2v 2 - 3) X 2 + (2a 2 v - §s<7i) (1 + v 2 )g^ X , (22) 

c = X V(l + i/ a ), (23) 

and v = icu/Q (so that unstable modes have Re(u) > 0). 
The discriminant is 

b 2 - 4ac = [(9 - Wis 2 ) (x + gSff 1 «rX 1 (l + f 2 )) 

-12(72a;V(l + v 2 )] [ X + isffiffl^l + ^ 2 )] . (24) 

<7|| does not appear in the dispersion relation and the lin- 
earized equations because k ■ B = 0. Thus the instability 
in an initially vertical magnetic field behaves identically in 
both the ambipolar diffusion (cri = 0, cm = 0) and Ohmic 
diffusion (cri = 0, cry = 02) limits. 



4.2 Ambipolar or Ohmic diffusion limit (<7i = 0) 

In this case (cf. Blaes & Balbus 1994) growth occurs up to 
a wave number 



k c v A /n = V3(l + X ~ 2 ) 



2\l/2 



(26) 



and the maximimum growth rate, again identified by requir- 
ing (J24) to vanish, is the positive real solution to 



12^ + 16x^o 



9x + 12 = . 



(27) 



For x ^ 1> vq is close to the ideal MHD value of |, declining 
substantially for x 1- For x. <tC 1, fo ~ §Xi an d fcoWA/f2 ~ 
v/fx- This is illustrated in the top panel of Fig. ^ which 
shows the growth rate as a function of wave number for x = 
00, 10, 1, and 0.1. There is not much change from the ideal 
limit until x < 10, when the growth rate and characteristic 
wave number decline rapidly as x is reduced. 



4.3 Hall limit (a 2 = 0) with a x B z > 

The departure from the ideal limit changes markedly when 
the Hall term is introduced. In the absence of ambipolar 
diffusion, with o\B z > 0, one finds growth for wavenumbers 
between and 



k c v A /Q= v / 3(l + 2 X 



-l\-l/2 



(28) 



and a maximum growth rate unchanged from the ideal limit, 
i.e. uo = § ! The corresponding wave number, however, de- 
creases as the geometric mean of Q and lo c as x "~ > 0: 



fc ^/fi = \/I(i + !x- 1 )" 1/2 



(29) 



This behaviour is demonstrated in the middle panel of Fig. 







4 RESULTS 

The dispersion relation ( po| ) is quadratic in k 2 , and for real v, 
c is positive. Thus there is at most one real, positive root to 
the dispersion relation. Typically the plots of v vs k resemble 
inverted quadratics, the growth rate increasing from zero 
at k = 0, rising to a maximum (denoted by uqQ.) at some 
wave number ko, and decreasing to zero at a critical wave 
number fc c . However, I show below that for a small region of 
parameter space, k c becomes infinite, and growth may occur 
over a wide range of wave numbers. 

4.1 Ideal MHD limit (x -> 00) 

Ideal MHD is recovered in the limit uj c 3> £1, in which case 
the dispersion relation reduces to 

(kv A /n) 4 + (2v 2 - 3)(kv A /tt) 2 + u 2 {l + v 2 ) = Q (25) 

As shown by Balbus & Hawley (1991), this yields growing 
modes for < kv A /Q, < v3, with the maximum growth rate 
i/o = | occuring when the discriminant of eq. ( j25|) vanishes, 
that is at kov A /Q — v^xf- 



4.4 Hall limit (a 2 = 0) with a x B z < 

If the other Hall-dominated limit is taken, the behaviour is 
again different (see lower panel of Fig. |l]). For \ > 2, the 
range of wave numbers for which growth occurs is finite, 
with 

k c v A /n = Vs(l-2 X ' 1 y 1/2 ■ (30) 

k c becomes infinite once x reaches 2, and for | < x < 2, 
all wave numbers grow. There is a fastest-growing mode for 
X > § , with i/ = § and 

fc <M/fi=yi(l-fx-T 1/2 (31) 

(cf. the x = 1-5 case in the lower panel of Fig. |l|). Again, this 
becomes infinite for x = § For i < x < |> the growth rate 
is a monotonically increasing function of k, asymptotically 
approaching 

^o = (x-§) 1/2 (2-x) 1/2 (32) 
This limit decreases to zero 

for X = 5 L cf - the X = 0.75 
and 0.51 curves in the lower panel of Fig. hh, at which point 
there are no longer any growing modes. 
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Figure 1. The growth rate of the Balbus-Hawley instability as a 
function of wavenumber for different choices of the Hall and Ped- 
ersen conductivities a\ and 02 (see text) . Upper: the ambipolar or 
Ohmic diffusion limits (cti = 0); middle: the Hall limit (<72 = 0) 
with a±B z > 0; and lower: the Hall limit for o\B z < 0. The 
curves in each panel are labelled by \, the ratio of the coupling 
frequency to the disc angular frequency: ideal MHD is recovered 
for \ — > 00, and the magnetic field is poorly coupled to the disc for 
X < 1. Note that the scale on the horizontal axis differs between 
the lower and upper two plots — the ideal MHD curve labelled 
"00" appears in each plot. 



4.5 The general case 

When o\B z > 0, there are growing modes with wave num- 
bers between and 



3x(x+ ttSQ-io-x 1 ) 
X 2 + f X s0 "! "! 1 + 1 



1/2 



(33) 



The maximum growth rate satisfies 
© 1999 RAS, MNRAS 000, §-§ 



Figure 2. As for Fig. [j], but for intermediate cases. Upper panel: 
o~\B z > and |cpl| = 02', middle panel: o\B z < and |ci| = 02; 
and lower panel: a\B z < with |<xi| = 61T2. 



12(T2<7 



V 9 



1 -1 



(l + ^o) =X- 



(34) 



This behavior is illustrated in the top panel of Fig. ^| which 
shows the wave number dependence of the growth rate when 
cri = <T2. As x is reduced, the maximum growth rate and 
associated wave number decrease. However, in the limit \ — > 
0, the growth rate remains finite, with 



The corresponding wave number, however, decreases: 

1/2 



3x 



(35) 



(36) 
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Figure 3. Maximum growth rate (in units of the Keplerian fre- 
quency H) as a function of degree of coupling between the mag- 
netic field and neutral gas, \, f° r different choices of the Hall 
and Pederson conductivities o\ and a^. Dotted curve: the am- 
bipolar and Ohmic diffusion limits (01 = 0); solid: the Hall limit 
(<T2 = 0) ; dashed: the case | o~i \ = 02 . The upper and lower mem- 
bers of the pairs of dashed or solid curves correspond to a\B z > 
or a\T3z < respectively. 



When a\B z < 0, growing modes exist only for \ > 
i|ai|/cr x . If I en I /<t ± > f , then eqs (||) and @ hold; oth- 
erwise there is a restricted range of \ for which all wave 
numbers grow, determined by the inequality 



5|oi 
4a ± 



5|oi 

4<T 1 



(37) 



The middle panel of Fig. shows the behaviour as \ i s re ~ 
duced for \o\\ = 02 . In this case the range of growing modes 
is always finite. The lower panel, for |ci| = 6(72 demonstrates 
the existence of a small range of \ around 0.75 over which 
all wave numbers grow. 

These results are summarised in Figures |^ and ^. Fig. ^ 
shows the maximum growth rate as a function of \ for differ- 
ent choices of <j\ and 02 ■ In the ambipolar and Ohmic diffu- 
sion limits (<ti = 0), the maximum growth rate decreases to 
zero as x ~~ * 0. When the Hall term is included, the growth 
rate drops to zero at x = |l cr i|/ (7 -L f° r v\B z < 0, or remains 
finite at j\ai\/(a± + a2) for aiB z > 0. This behaviour is dis- 
played more fully in Fig. ^| which shows contours of fixed 
growth rate in a polar plot. The ambipolar and Ohmic dif- 
fusion limits apply along the locus o\ = running vertically 
through the origin. It is clear that the growth rate in the 
weak coupling limit (x Ss 1) is sensitive to the form of the 
conductivity tensor, and that the behaviour in the two dif- 
fusion limits is not generally representative. 



5 A TWO-FLUID MODEL FOR THE HALL 
LIMIT 

It is surprising at first glance that the Hall conductivity 
allows the Balbus-Hawley instability to grow at a rate of or- 
der fi in the formal limit of zero coupling. Here I address the 
physical explanation for this result, for simplicity focussing 
on the "pure" Hall case {02 = 0). In the ambipolar diffu- 
sion limit, intuition is aided by a model two-fluid system: an 



D 




-2-10 1 2 

sign(B z )%o 1 /o L 

Figure 4. A polar plot of the maximum growth rate of the 
Balbus-Hawley instability. The Cartesian cordinatcs of a point 
on the plane corresponds to (sign(B z )oi , 02) suitably normalised. 
The polar radius is \, the parameter controlling the degree of 
magnetic coupling in the disc. The contours are labelled by the 
maximum growth rate in units of (i~ 1 . 



ionized (but electrically-neutral) fluid into which the mag- 
netic field is frozen; and a neutral fluid which is indirectly 
coupled to the magnetic field via collisions of the charged 
speicies with the neutrals. There is an analogous model in 
the Hall limit |<ti| 3> 02- 

In this model, the ionized component of the fluid can 
be regarded as consisting of two species, one of which does 
not suffer significant collisions with the neutrals, the other 
being so strongly tied to the neutrals by collisions that its 
drift velocity is negligible. For convenience, I shall refer to 
the former as "electrons" (subscript e) the latter as "ions" 
(subscript i); although the roles of these species could be 
played by other species, e.g. ions and grains, under suitable 
conditions. 

The Hall parameters of the two species satisfy \/3i\ <C 
1 <C |/3 e |. Because the ions are strongly tied to the neutrals 
by collisions and are responsible for transmitting the elec- 
tromagnetic stresses to the fluid, one may regard the ions 
and neutrals as a single fluid with the charge density of the 
ions and the inertia (and thermal pressure, when relevant) 
of the neutrals. This fluid is permeated by the electrons, into 
which the magnetic field is frozen. 

In this limit the components of the conductivity per- 
pendicular to the magnetic field are 



and 



<T-2 



ecriiZi 
B 



(38) 



(39) 



(where I have used riiZi = —n e Z e in eq. |38|). The current is 
dominated by the electrons, which drift perpendicular to the 
magnetic and electric field so that the net Lorentz force on 
them is zero - collisions with the neutrals play a negligible 
role in their equation of motion. The drift of the ions with 
respect to the neutrals, however, is strongly inhibited by 
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collisions, and is determined by the balance between the 
electric field and the drag associated with the neutrals. Thus 

J « n e Z e ev e « a ± B X E± . (40) 

The electromagnetic stresses on the fluid are communicated 
largely by the collisions with the ions: 

« aiE ± « jiPPiVi (41) 

c 

Because the ion-neutral drift speed Vi is small, the associ- 
ated dissipation is small: J ■ E' ± w o^E™, tending to zero 
in the limit 01 — * 0. With this model fluid, one can use the 
linearized equations (p^)-(p^) to understand the growth of 
the instability in the Hall limit as \ ~~ > 0. 

First consider the wave modes supported by the fluid 
in the absence of rotation. At small wavenumbers (long 
wavelengths) the fluid supports the Alfven waves of ideal 
MHD, while at wavenumbers far in excess of lo c /va, there 
are two circularly-polarized modes (Wardle & Ng 1999). In 
one, the electrons and magnetic field are relatively unper- 
turbed, and the ion-neutral fluid executes gyrations about 
the field direction at the frequency uj ~ U) c . In the other, 
the ion-neutral fluid is static and the electrons and mag- 
netic field twist about the unperturbed field direction with 
the opposite sense of rotation. The frequency of this mode 
scales quadratically with wavenumber, as the diffusive term 
V X E' dominates V X (v X B) in the induction equation: 

{kv A ) 2 /uj c . The sense of rotation of either mode is de- 
termined by the sign of cr\B z . 

Finally, consider the effect of rotation, which attempts 
to enforce epicyclic motion of the ions and neutrals with 
frequency Q through the Coriolis and angular terms in the 
momentum equation. For \ ^ 1j ^ 2> w o, so the rota- 
tion can only couple effectively to the second mode for wave 
numbers kvA ~ \fuj c Sl. The generation of toroidal field from 
poloidal field by the shear in the disc (see eq. ^J) acts as a 
forcing term for this circularly polarized mode, which grows 
or decays depending on whether the shear assists or opposes 
the rotation of the magnetic field, i.e. depending on the sign 
of o\B z . 



6 DISCUSSION 

The Hall conductivity modifies the growth of the Balbus- 
Hawley instability in weakly ionized discs when the mag- 
netic field is not well-coupled to the neutral matter. The 
effect depends on the relative signs of o\ and the initially 
vertical field B z . For o\B z > the growth rate remains of 
order the Keplerian frequency Q as the coupling is reduced, 
and the wavelength of the fastest growing mode scales as 
X 1 / 2 . On the other hand, if a±B z < 0, the mode is damped 
for x 5s \- The dependence on the initial field being parallel 
or antiparallel to the disc rotation axis arises ultimately from 
the intrinsic handedness introduced into the medium by the 
difference in the drifts of positive and negative species. In the 
limit SI — > 0, the unstable mode corresponds to a standing, 
circularly-polarized Alfven wave; left and right-circularly po- 
larized Alfven waves propagate differently according to the 
sign of g x B z (Wardle & Ng 1999). 

This is to be contrasted with the ambipolar and Ohmic 
diffusion limits (Blaes & Balbus 1994; Jin 1996), where 



the growth rate is independent of the sign of B z , declin- 
ing steadily while the wavelength increases linearly as the 
coupling parameter is reduced. Blaes and Balbus found that 
the growth rate increases again for x ~> Pi/ Pi when the cou- 
pling becomes so poor that the momentum-exchange time 
scale for the ions with the neutrals^, (jip) -1 is longer than 
the Keplerian time scale f2 _1 . Then the instability develops 
in the ionized component and magnetic field without affect- 
ing the neutrals. This behaviour does not occur here because 
I neglected the inertia of the charged species (see Section ) 
- the formulation in Section ^| is valid only for frequencies 
well below jip, and thus implicitly assumes that Q <C ytp, 
i.e. x ^ Pi/ ' P- The ionization fraction in protostellar discs 
is so low that this approximation is very good indeed, and 
if the coupling were this poor any dynamically important 
magnetic field would remove itself (and the ions) from the 
disc on the Alfven crossing time in the ionized fluid. The ion- 
ized species cannot, therefore, decouple from the neutrals in 
any practical sense. 

For a favourable alignment between the magnetic field 
and rotation axis (i.e. o\B z > 0), the growth rate in the 
weak-coupling limit remains within a factor of a few of Q. 
The practical limit on the coupling is determined by the 
wavelength as the coupling gets weak, as this must be less 
than the scale height h of the disc. If c s is the isothermal 
sound speed, so that h = c s /f2, the wavelength of the fastest- 
growing mode becomes of order h when the coupling pa- 
rameter x reaches v\/c? a . The corresponding limit in the 
ambipolar diffusion case is va/c s which would produce a 
growth rate of order (va/c 3 )SI. 

A study of the nonlinear development of the instabil- 
ity requires a numerical simulation. For now, I note that 
the strength of the magnetic field affects the conductivity 
through the Hall parameters, with a larger field increasing 
the degree of flux- freezing in the ionized component. Thus 
there is the possibility of the instability beginning in the Hall 
regime, with the field strength growing to the point that the 
fluid enters the ambipolar diffusion regime. 

These results are applicable to protostellar discs, where 
collisions are sufficient to decouple ions and grains from the 
magnetic field (Wardle & Konigl 1993). Indeed, the Hall 
conductivity is of the same order as the ambipolar diffu- 
sion term over a broad range of parameters (Wardle & Ng 
1999). Further the coupling parameter is thought to be small 
near the midplane in protostellar discs, particularly at radii 
around one AU where the temperature is low enough that 
thermal ionization is not yet effective, and the surface den- 
sity is high enough to shield the material from cosmic rays 
(Hayashi 1981; Gammie 1996; Wardle 1997). The magni- 
tude and relative size of the components of the conductivity 
tensor vary strongly with height in the disc as the neutral 
density decreases near the disc surface, and the ionization 
level increases because of exposure to cosmic rays (Gammie 
1996; Wardle 1997) and X-rays from the central star (Glass- 
gold et al 1997). In particular, although the instability is not 
very efficient in transporting angular momentum when the 
coupling is weak, the column density within this portion of 



* As distinct from the momentum exchange time scale for the 
neutrals with the ions, (~ipi)~ 1 that appears in the definition of 
X in the ambipolar diffusion limit. 
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the disc is relatively large, so that the net angular momen- 
tum transport can be comparable to that in the layers closer 
to the surface. 

The dependence on the sign of the initial field raises 
the intriguing possibility that the accretion process, and the 
formation of magnetically-driven disc winds are dependent 
on the sign of the magnetic field. This issue is linked to the 
collapse of cloud cores to form protostars and discs, which 
presumably brings in some fraction of the interstellar field 
to be the initial disc field. While the calculations of core 
collapse to date are instructive (e.g. Fiedler & Mouschovias 
1993; Basu 1997; Contopoulos, Ciolek, & Konigl 1998; Li 
1998; Ciolek & Konigl 1998), none include the Hall effects, 
which are also important during this phase (Wardle & Ng 
1999). Angular momentum transport during collapse will 
also depend on whether the field is parallel or antiparallel to 
the core's angular momentum vector, and one might imag- 
ine that this leads to a preferential orientation within discs, 
but this is unclear as the sign of a\ may change between the 
collapse and disc phases as grains, then ions become suc- 
cessively decoupled from the magnetic field with increasing 
neutral density. 



7 SUMMARY 

This paper examined the Balbus Hawley instability of a ver- 
tical magnetic field in a weakly-ionized, thin Keplerian disc. 
The vertical stratification was ignored, and radial and az- 
imuthal variations were neglected. The resulting linearized 
equations are valid for wavelengths small compared to the 
disc scale height. I examined the dispersion relation. The 
conductivity tensor contributes two parameters: (i) a cou- 
pling parameter the ratio of the frequency at which the 
finite conductivity becomes important to the Keplerian an- 
gular frequency; and (ii) the ratio sign(B z )ai/crj_, which re- 
flects the contributions of the Hall and Pedersen conductiv- 
ities, oi, (72 respectively, to the total conductivity perpen- 
dicular to the magnetic field a± = [p\ + u 2 ) 1 ^ 2 ■ 
I found the following: 

(i) The Hall conductivity substantively modifies the be- 
haviour of the instability once finite conductivity effects be- 
come important (x & 1). If the rotation axis of the disc 
defines the z direction, the departure from the ambipolar 
diffusion limit depends on the sign of o\B z . 

(ii) For <J\B Z > 0, growth occurs over a finite range of 
wavenumbers. For poor coupling (x <§C 1), the maximum 
growth rate remains finite, tending to a value §|o"i|/(<jx + 
02) while the associated wave number scales as ^f\. This 
should be contrasted with the ambipolar diffusion case 
(Blaes & Balbus 1994) for which the growth rate and wave 
number scale linearly in \- 

(iii) For a±B z < 0, on the other hand, there is no instabil- 
ity when the coupling is poor: growth requires \ > kWl \/&± ■ 
There is generally a limited range of wave numbers, over 
which growth occurs. However, if |ai|/oj_ < |, there is a 
limited range of \ around the value \a-i_\/a± for which all 
wave numbers grow. 

(iv) Although the instability has a formal growth rate of 
order il in the limit of zero coupling when aiB z > 0, the 
wavelength is limited by the disc scale height, and there- 
fore growth will only occur for \ ii v a/c 2 - Note that this 



constraint is far less severe than in the ambipolar diffusion 
limit, which requires \ iS va/c s - 

Finally, I stress that these results indicate the qualitative 
changes to dynamical behaviour to be expected when am- 
bipolar diffusion breaks down during the collapse of cloud 
cores to form protostars and during the subsequent evolu- 
tion of their attendant discs. 

I thank Charles Gammie and the referee for comments 
on the manuscript. The Special Research Centre for The- 
oretical Astrophysics is funded under the Special Research 
Centres programme of the Australian Research Council. 
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